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We study the pairwise quantum discord (QD) for a symmetric multi-qubit system in different 
types of noisy channels, such as phase-flip, amplitude damping, phase-damping, and depolarizing 
channels. Using the QD and geometric measure of quantum discord (GMQD) to quantify quantum 
correlations, some analytical or numerical results are presented. The results show that, the dynamics 
of the pairwise QD is related to the number of spin particles N as well as initial parameter 6 of 
the one-axis twisting collective state. With the number of spin particles N increasing, the amount 
of the pairwise QD increases. However, when the amount of the pairwise QD arrives at a stable 
maximal value, the pairwise QD is independence of the number of spin particles N increasing. The 
behavior of the pairwise QD is symmetrical during a period 0 < 0 < 27r. Moreover, we compare 
the pairwise QD dynamics with the GMQD for a symmetric multi-qubit system in different types 
of noisy channels. 
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I. INTRODUCTION 

As we all know, quantum correlations, such as quan¬ 
tum entanglement, have been proposed as the key re¬ 
source present in certain quantum communication tasks 
and quantum computational models [1-2]. Many different 
quantum correlation measures have been proposed to de¬ 
tect the nonclassical correlations beyond entanglement. 
Among them, quantum discord (QD) which has initially 
been introduced by Olliver and Zurek [3], captures more 
general quantum correlations, and later the studies of 
QD have been vigorously investigated both theoretically 
[4-11] and experimentally [12,13]. However, evaluation 
of QD requires a potentially complex optimization pro¬ 
cedure in general and analytical results have been obtain 
only in few restricted cases of 2 0 2 and 2 0 A^ systems 
[14-22]. In order to overcome the optimization proce¬ 
dure of QD, Dakic et al. [23] introduced another measure 
quantum correlations named geometric quantum discord 
(GMQD)for a more general bipartite quantum system. 

On the other hand, decoherence is an inevitable phe¬ 
nomena owe to quantum system inevitable interacting 
with its surrounding environment, which gives rise to 
degradation of quantum coherence. In other word, deco¬ 
herence is a major hindrance to practice quantum infor¬ 
mation processing. Hence, it is very interesting and nec¬ 
essary to study the effect of decoherence on the evalua¬ 
tion of quantum correlations. Recently the researches on 
dynamics of pairwise entanglement have been made great 
progress in various fields [24-32]. For example, Wang 
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et al. showed spin squeezing implies pairwise entangle¬ 
ment for arbitrary symmetric multi-qubit states [28,29]. 
Yin et al. investigated the pairwise quantum correlations 
of symmetric multi-qubit states by using the geometric 
measure of quantum discord [30]. Xi et al. investigated 
pairwise quantum correlations for superpositions of Dicke 
states [32]. But pairwise correlations, as we know, quan¬ 
tified by QD and its geometric measure, in a symmet¬ 
ric multi-qubit system under decoherence noisy channels 
seem to have been seldom exploited before. It is of inter¬ 
est and importance to investigate how decoherence affects 
the pairwise correlations under decoherence noisy chan¬ 
nels. Motivated from the recent study on the dynamics 
of QD and GMQD under the influence of external en¬ 
vironments for different multipartite states [33], such as 
Werner-GHZ type three-qubit and six-qubit states, we 
here devote to examining the pairwise quantum corre¬ 
lation properties in terms of QD for pairs of particles 
extracted from a symmetric state. The two-particle den¬ 
sity matrix is expressed in terms of expectation values of 
collective spin operators S for the large system [28-30,34- 
40]. 

Our attention is mainly focused on what happens to 
the pairwise QD for a symmetric multi-qubit system in 
four different types of noisy channels, such as phase- 
flip, amplitude damping, phase-damping, and depolar¬ 
izing channels. We extend the results obtained in Ref. 
[40], where using the spin squeezing to witness multipar¬ 
tite correlations under decoherence is investigated. Using 
QD and GMQD to quantify quantum correlations, some 
analytical or numerical results are presented. The results 
show that, the dynamics of the pairwise QD is related to 
the number of spin particles N as well as initial param¬ 
eter 9 of the one-axis twisting collective state. With the 
number of spin particles N increasing, the amount of the 
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pairwise QD increases. However, when the amount of the 
pairwise QD arrives at a stable maximal value, it is in¬ 
dependence of the number of spin particles N increasing. 
The behavior of the pairwise QD is symmetrical during 
a period 27r. Moreover, we compare the pairwise QD 
dynamics with the geometric measure of GMQD for a 
symmetric multi-qubit system in different types of noisy 
channels. 

The layout is as follows: In Sec. II, we briefly review 
the definition of QD and GMQD, and give the expres¬ 
sion of QD and GMQD. We illustrate initial states and 
noise channels in Sec. III. We devote to examining their 
pairwise quantum correlation properties in terms of QD 
for pairs of particles extracted from a symmetric state of 
multi-qubit systems in different types of noisy channels 
in Sec. IV. Finally, we give the conclusion in Sec. V. 

II. QD AND GMQD EXPRESSION OE TWO 
QUBITS 

Before starting to discuss the dynamics of quantum 
correlations, we review the correlation measures used in 
our investigation, namely, QD and GMQD. Let us start 
with the definition of QD for any bipartite quantum X- 
state pAB 

/ Pii 0 0 Pi4 

_ 0 P22 P23 0 

0 P32 P33 0 

V P41 0 0 P44 

then, the classical correlation C{pab) is defined as 

C{pab) = S{pa) - min S{pa\b), (2) 

mi} 

where S{p) = —tr{p\og 2 p) is the von Neumann entropy. 
Moreover, the minimum is taken over the set of posi¬ 
tive operator valued measurement {nf} on subsystem 
B, S{pa\b) is the conditional entropy for the subsystem 
A. QD can be simply obtained by subtracting C{pab) 
from the total amount of correlation [3] 

QD{pab) = Hpa-.b) — C{pab), (3) 

where the total correlation is quantified by the quantum 
mutual information I{pa-.b) 

Hpa-.b) = S{pa) + S{pb) + T,i^ie^\og2ei, ( 4 ) 

where pa(pb) is the reduced matrix of pab by tracing 
out B{A), and 

S{pa) = —(pii + P 22 ) log 2 (pii + P 22 ) 

- (P33 + P 44 ) l0g2(P33 + P 44 ), (5) 

S{pb) = —{pil + P 33 ) log 2 (pll + P 33 ) 

- {P22 + Pn) log2(p22 + P44)- (6) 


The eigenvalues of the density matrix pab in Eq. (II) 
are given by 


El — 5[(pil + P44) + \/(pil ~ P 44 )^ + 4 |pi 4 p], 

£2 = \ + P 44 ) - \ /(pll -P44)^ +4|pi4p ], 

£3 = ^[(^22 + P 33 ) + \ /{P22 — P33)^ + 4|p23p ]; 

£4 = \ [{P22 + P 33 ) — a/ {P22 — P 33 Y + 4|p23p]- 

Since the calculation of classical correlation involves 
a potentially complex optimization process, there exists 
no general analytical expression of discord even for the 
simplest case of two-qubit state. It is difficult to calcu¬ 
late the QD by analytical solutions, while for a bipartite 
quantum X-state described by the density matrix pab, 
we can derive an expression of the QD [22] 

QD{pab) = min((5i, ( 52 ), (8) 

where 

Qi = H{pii + P33) -I- log2 £i -I- Dj, 

_ jj-j- l-t-Vfl — 2(p33-eP44)]^+4(|pi4|-|-|p23l)^ ^ 

D 2 = -T,iPii log 2 Pii - H{pii + P 33 ), 

H{x) = -a;log 2 X - (1 - x) log 2 (l - x). 

Besides, GMQD for any bipartite quantum X-state 
PAB bas been proposed by Dakic et al. [23]. A bipar¬ 
tite quantum X-state pab writes in Bloch representation 
as 

I ^ 

PAB = -^[o-Q G CTo + G CTo -I- 2/iCTo G CTi) 

i 
3 

+ ^ Rij(7i<Si(Jj], (9) 

i,i=l 

where Xi = Tr{pai G cto), yi = Tr{pao G <7^) are com¬ 
ponents of local Bloch vectors and Rij = Tr{pai G (Jj) 
are components of the correlation tensor. The operators 
ai{i = 1,2,3) stand for the three Pauli matrices and (Tq 
is the identity matrix. The expression of the GMQD is 
given by 

GMQD{pab) = i(||xll" + \\Rf - k^ax), (10) 

where x = {xi,X 2 ,X 3 )'^, R is the matrix with elements 
Rij , and kmax is the largest eigenvalues of matrix defined 
by 

K = xx^ + RR'^. (11) 

Note that the function GMQD (pab) reaches its maxi¬ 
mal value of 1/2 , the normalized GMQD takes the form 
2GMQD{pab) in this paper. In the following, we use 
Eq. (8) and Eq. (10) to discuss the pairwise quantum 
correlation dynamics under different types of noisy chan¬ 
nels. 
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III. INITIAL STATES AND NOISE MODEL 

We consider a system of N exchange symmetry spin- 
1/2 particles with the ground state |0) and excited |1). 
This system’s properties can be described by the col¬ 
lective operators Sa = fo’’ 

a = x,y,z. The one-axis twisting Hamiltonian reads[40] 

H = xSl (12) 

where x is a nonlinear coupling constant for all pair¬ 

wise particles interaction. We choose the initial symme¬ 
try state, which has been prepared in the product state 
\0)n = |00...0), based on its dynamic evolution, the one- 
axis twisting collective state with even parity at time t is 
obtained as 

IV'(t)) =exp(-i6»5'^)|0)jv, (13) 

where 9 = xt is the one-axis twisting angle. Following 
the previous discussion, the states with exchange symme¬ 
try and parity ensures that its two-qubit reduced state 
can be extracted randomly from this state, in the basis 
of |00), |01), |10), |11), whose two-quibt reduced density 
matrix can be written form Eq. (1). The elements of 
the two-quibt reduced density matrix can be represented 
by the local expectation values for the one-axis twisting 
state [28-30,34-40] 

Pii = ^(1 + 2 (cti2) -I- {aiza2z)), 


P22 = P33 = ^(1 - {<Jlz(J2z)), 


P44 = ^(1 - 2{aiz) + {(7iza2z)), 

P23 = P 32 = (o'H-0'2-), 

Pl4 = P*41 = (cTl-CTa-). (14) 

where the local expectation values for the one-axis twist¬ 
ing state were given by 

(o-iz) = -cos^-^(^). 


(cri_cr2-) =-i(l - cos^ ^0) - ^sin(^)cos^ ^(^). 

(15) 

Following, we will focus on decoherence of the one-axis 
twisting state when each spin particle is independently 
coupled to local noisy environment. Decoherence is a ma¬ 
jor hindrance to practice quantum information process¬ 
ing, it is important to discuss how decoherence affect the 
quantum correlation. Here, we investigate what happens 
to the pairwise QD for a symmetric multi-qubit system 
in four different types of noisy channels, such as phase- 
flip, amplitude damping, phase-damping, and depolariz¬ 
ing channels. The dynamics of qubits interacting inde¬ 
pendently with individual environments is described by 
the solutions of the appropriate Born-Markov-Lindblad 
equations, which can also be obtained by so called the 
Kraus operator approach. Given the initial state p, the 
state under noise is given by 

N N 

<p)= E (16) 

E1,E2,..EN i=l i=l 

where Kei denotes the Kraus operator for the z—th par¬ 
ticle. Since the system has exchange symmetry, and the 
decoherence channels act independently on each parti¬ 
cle. Therefore, all the elements of the two-quibt reduced 
density matrix are determined by some correlation func¬ 
tions and expectations, and any expectation value of the 
operator A can be calculated as (A) = Tr[Ae{p)] = 
Tr[si{p)A]. 

A. Phase flip channel 

The Kraus operators for single qubit in phase-flip chan¬ 
nel are given by 

= (U) 

E. = (f _[/;). (18) 

where p = exp(—qt). Since the system has exchange sym¬ 
metry, and the decoherence channels act independently 
on each particle. Therefore, substituting Eqs. (17) and 
(18) into Eq. (16), and combining Eq. (14) and (15), the 
elements of the density matrix can be represented by the 
local expectation values 

Pii(f) = ^(1 + 2(o-iz) -f {aiza2z)), 


{<Jlz'J2z) ^{l+COS^ "^9), P 22 (f) = P 33 = ^(1 - (ct 1 zC^ 2 z)), 

(CTi + Cr2-) = -^(1 - COS^“^ 9), P44{t) = ^(1 “ 2((Tlz) -f (cri^Cr2z)), 



4 


P23it) = PaaW = (1 - 2 p)2(cti+CT2-), 

Pl4it) = pli{t) = (1 - 2p)2(cri_CT2-). (19) 

B. Amplitude damping channel 

Amplitude damping channel which is used to charac¬ 
terize spontaneous emission describes the energy dissipa¬ 
tion from a quantum system. The Kraus operators for a 
single qubit are given by 




( 20 ) 


where p = exp(—yt). According to previous analysis, the 
elements of the density matrix can be represented by the 
local expectation values 

Pll{t) = + 2(o-lz) -b {(J1 z(T2z)), 


P 22 {t) = P33 = ^(1 - (crizCT 2 z)), 


Piiit) = -(1 - 2{aiz) + {(Tiz(J2z)), 

P23{t) = P*32it) =/(0-l-rCr2-), 
Pl4it) = P*4i{t) =/(cri_cr2-)- 


(26) 




0 0 
0; ’ 


( 21 ) 


where p = exp(—yt). According to previous analysis, the 
elements of the density matrix can be represented by the 
local expectation values 

Pii{t) = i[l -b 2{p{aiz) -bp - 1) -b ip^{<7iza2z) 

- 2(1 -p)p(cti^) -b (1 -p)^)]. 


P22(t) = P33 = -[1 - {p^{(Jiza2z) - 2{1 - p)p{aiz) 

+ (l-pf)], 

P4-i{t) = ^[1 “ 2 ( p ( ctu ) + p- 1 ) + {p^{aiz(T2z) 

- 2(1 -p)p(cri^) + (1 -p)^)], 

P23{t) = P32{t) =p{ai + a2-), 

Pl4{t) = P4i(t) = p(cti_(T2-)- (22) 

C. Phase-damping channel 

Phase-damping channel describes a quantum noise 
with loss of quantum phase information without loss of 
energy. The Kraus operators for a single qubit are given 

by 


En — 


E, = 


E 2 — 


Vp 0 
0 Vp 


0 

0 0 


0 0 

0 


D. Depolarizing channel 

Depolarizing channel is another important type of 
quantum noise, it describes the process in which the den¬ 
sity matrix is dynamically replaced by the state 1/2. I 
denoting identity matrix of a qubit. The Kraus operators 
for a single qubit are given by 


An — A/I — 


3p / I 0 

0 1 




E 2 — 


IP f 0 -i 

4 V i 0 


E3 = A - 


P f I 0 


4 ( 0-1 


(27) 

(28) 

(29) 

(30) 


where p = exp(—y<). According to previous analysis, the 
elements of the density matrix can be represented by the 
local expectation values 

Pllit) “ — p){<^lz) + |) + (1 — p)^(CTizCr 2 z) 

+ (1 -p)p(ctiz) + ^], 

P44{t) = ^[1 - 2((1 -p)(criz) + |) + (1 -p)^(CTlzCr2z) 

+ (1 -p)p(ctiz) + ^], 


(23) 

P22{t) — P 33 — ^[1 - ((1 

(24) 



P23{t) = Phit) 

(25) 

Pl4{t) = P*4i{t) 


(31) 
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IV. DISCUSSION 

In this section, we devote to examining their pairwise 
quantum correlation properties in terms of QD for pairs 
of particles extracted from a symmetric state of multi¬ 
qubit systems in different types of noisy channels, such 
as phase flip, amplitude damping, Phase-damping, and 
depolarizing channel, the two-particle density matrix is 
expressed in terms of expectation values of collective spin 
operators S for the large system. Combining Eq. (18)- 
Eq. (20), we can compute QD. 

In Figure 1, we plot the pairwise QD of the one-axis 
twisting collective spin state given by Eq. (11) as a func¬ 
tions of yt and N for fixed 9 — O.Itt in different types of 
noisy channels. It can be seen that in phase-flip chan¬ 
nel, the dynamics of pairwise QD decays quickly to zero 
and revives increasingly to its initial values. This result 
agrees with previous study of similar nature for GHZ 
states [33]. The revival of pairwise QD dynamics reaches 
to its initial values and keeps this stable values in phase- 
flip decoherence channel. This indicates the pairwise QD 
is not effected by phase-flip decoherence channel after 
The revival of pairwise QD dynamics reaches to its ini¬ 
tial values. Besides, with the number of spin particles 
N increasing, the amount of the pairwise QD increases. 
However, when the amount of the pairwise QD arrives at 
a stable maximal value, it is independence of the number 
of spin particles N increasing. Compared to phase-flip 
channel, the pairwise QD dynamics decay in a monotonic 
fashion in amplitude damping, phase-damping, and de¬ 
polarizing channel. This indicates that the amplitude 
damping, phase-damping, and depolarizing channel are 
heavily influence the pairwise QD. What is more, the 
pairwise QD dynamics decays more slowly in amplitude 
damping channel than those in Phase-damping, and de¬ 
polarizing channel. 

Next, we investigate the influence of initial angle 6 of 
the one-axis twisting collective state on the dynamics of 
the pairwise QD under different types of noisy channels. 
We plot the pairwise QD of the one-axis twisting collec¬ 
tive spin state given by Eq. (11) as a functions of yt and 
6 for fixed N = 12 in different types of noisy channels, 
as shown in Figure 2. It is clearly seen that, in phase- 
flip channel, the pairwise QD quickly decays to zero and 
revives. When the amount of the pairwise QD arrives 
at its initial value, it is not effected by phase-flip noisy 
channel. This indicates that QD for pairs of particles 
extracted from a symmetric state of multi-qubit systems 
can persist at a steady state. Besides, the behavior of the 
pairwise QD is symmetrical during a period 2tt. However, 
compared to phase-flip channel, the pairwise QD dynam¬ 
ics decay in a monotonic fashion in amplitude damping. 
Phase-damping, and depolarizing channel. 

Finally, we compare the pairwise QD dynamics with 
the geometric measure of quantum discord (GMQD) for 
a symmetric multi-qubit system in different types of noisy 
channels. Figure 3, presents dynamics of QD and GMQD 
for the one-axis twisting collective spin state given by 


Phase - flip 



Amplitude - damping 



Depolarizing - damping 



FIG. 1. (Color online)Quantum discord (QD) of the one-axis 
twisting collective spin state given by Eq. (11) as a functions 
of yt and N for fixed 0 — O.Itt in different types of noisy 
channels. 
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FIG. 3. (Color online) Dynamics of quantum discord (QD) 
and geometric measure of quantum discord (GMQD) for the 
one-axis twisting collective spin state given by Eq. (11) with 
9 — O.Itt, a = 12 in different types of noisy channels. 


yt 


Depolarizing-damping 


QD 


V'^^GMQD 
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FIG. 2. (Color online)Quantum discord (QD) of the one-axis 
twisting collective spin state given by Eq. (11) as a func¬ 
tions of 'yt and 6 for fixed A = 12 in different types of noisy 
channels. 
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Eq. (11) with 0 = O.Itt, TV = 12 in different types of 
noisy channels. It is also observed that the pairwise 
QD has subtle differences with GMQD. Quantum cor¬ 
relations quantified by QD are always larger than those 
quantified by GMQD in decoherence noisy channels ex¬ 
cept for depolarizing noisy channel. This means the QD 
is more robust than GMQD under the influence of deco¬ 
herence noisy channels, and depolarizing noisy channel 
heavily influences the pairwise quantum correlations dy¬ 
namics. 


V. CONCLUSION 

We study the pairwise QD for a symmetric multi-qubit 
system in different types of noisy channels, such as phase- 
flip, amplitude damping, phase-damping, and depolariz¬ 
ing channels. Using QD and GMQD to quantify quan¬ 
tum correlations, some analytical or numerical results are 
presented. The results show that, the dynamics of the 


pairwise QD is related to the number of spin particles 
N as well as initial parameter 9 of the one-axis twist¬ 
ing collective state. With the number of spin particles 
N increasing, the amount of the pairwise QD increases. 
However, when the amount of the pairwise QD arrives at 
a stable maximal value, it is independence of the number 
of spin particles N increasing. The behavior of the pair¬ 
wise QD is symmetrical during a period 27r. Moreover, we 
compare the pairwise QD dynamics with the geometric 
measure of GMQD for a symmetric multi-qubit system 
in different types of noisy channels. It is perhaps helpful 
for the quantum information processing. 
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